Disorder Induced Quantum Phase Transition in Random-Exchange Spin-1/2 Chains 
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We investigate the effect of quenched bond-disorder on the anisotropic spin-1/2 (XXZ) chain 
as a model for disorder induced quantum phase transitions. We find non-universal behavior of 
the average correlation functions for weak disorder, followed by a quantum phase transition into 
a strongly disordered phase with only short-range xy-correlations. We find no evidence for the 
universal strong-disorder fixed point predicted by the real-space renormalization group, suggesting 
a qualitatively different view of the relationship between quantum fluctuations and disorder. 
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The existence and nature of quantum phase transi- 
tions (QPTs) |lj H, U has in recent years emerged as 
one of the most interesting aspects of low-dimensional 
quantum systems. QPTs arise from the subtle inter- 
play between short-range interactions on one hand and 
quantum fluctuations on the other Q. Since the lat- 
ter are particularly strong in one dimension, quantum 
spin chains have emerged as a generic model to inves- 
tigate QPTs H g, @, §. The additional presence of 
disorder has profound effects on the properties of low 
dimensional systems |], [Xo| as it competes with the sub- 
tle effects of quantum fluctuations. Its effect on QPTs 
has been the subject of recent intense and controversial 
discussion @,[§p[i|[D§0j^|l§. 

In one-dimension the strong-disorder renormalization 
(SDRG) (l7[ [l8) group offers potentially exact results 
for a variety of models. Of particular interest is the 
prediction of a universal infinite randomness fixed point 
(IRFP) for disordered spin chains. In many systems 
SDRG studies suggest the relevance a random- singlet 
(RS) phase|n], as the ground state for fairly general 
disorder. In this RS phase the average spin correlations 
are predicted to obey a universal isotropic power-law de- 
cay, KSfSj 1 )] ~ \i — j\~ 2 (a — x,z), where the over- 
bar denotes a configurational average over many random 
chains (replicas). The Luttinger (or spin- liquid) contin- 
uum of critical ground states of the ordered chain is thus 
predicted to collapse to a single point. Numerical re- 
sults consistent with the RS picture were reported for 
relatively short (N < 18) XXZ chains Q and also for 
long XX chains (A = 0) [|l3| with couplings uniformly 
distributed in [0,1]. Recently, some studies suggest the 
relevance of such a fixed point for the q-state quantum 
clock model and the quantum Ashkin Teller model [B|, 
while others dispute its existence, at least for S > 1/2 [3|. 

Recent numerical advances, in particular the develop- 
ment of the density matrix renormalization group 
now offer a framework to investigate the relevance of the 
IRFP to realistic one-dimensional spin systems p3j. In 
this Letter we investigate the influence of exchange dis- 
order on the anisotropic spin 1/2 Heisenberg chain (XXZ 



model), one of the best-known model systems for QPTs 
in one dimension. We find a qualitatively different sce- 
nario for the interplay of quantum fluctuations and dis- 
order. Our results indicate that the spin correlations do 
not obey the universal parameter independent decay law 
suggested by the RS picture. Instead we find a disorder- 
induced QPT for finite disorder strength, whose nature 
can be illustrated by an exactly solvable model. 

In this paper we present results of a density matrix 
renormalization group study of XXZ chains with ran- 
domness in the transverse nearest-neighbor coupling, 

N—l 

H = J £ [\i(S*Sf +1 + SfSf +1 ) + ASfSfa] , (1) 

(J > 0) where the anisotropy parameter A > controls 
the relative strength of the quantum fluctuations. In the 
homogeneous system (Ai = 1) the ground state of (|l|) 
shows long-range order for A > 1 (Ising regime) , whereas 
for A < 1 (critical regime) the spin correlations decay 



algebraically to zero as \{S ^Sj )\ 
universal decay exponents |23| 



Va with non- 



Vx = Vz 
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The introduction of quenched randomness brings 
about subtle changes in the ground-state. For reasons 
of better numerical control of replica averages we used a 
bounded probability density p(X) — ^r&(W — |A — 1|). 

We investigated chains of length up to L — 80 with a 
finite-size DMRG algorithm and ensured that the ground 
state and correlations (SfS") could be determined with 
consistent accuracy for arbitrary choices of A^. We noted 
that for long chains the standard DMRG procedure tends 
to spontaneously break the local S z reversal symmetry. 
We ensured that the resulting data for the correlation 
functions were nevertheless correct by comparing with 
explicitly symmetry-adapted high-accuracy DMRG cal- 
culations for a number of randomly selected replicas Xi. 
We verified the accuracy of the correlation functions and 
energies by comparing with exact data for XX (A = 0) 
chains which can be mapped to non-interacting lattice 
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FIG. 1: Decay exponents r\ x (circles) and rj z (squares) for 
A < 1 obtained with the DMRG for L=80 in comparison with 
analytic results for the thermodynamic limit. OBC and PBC 
denote open and periodic boundary conditions, respectively. 
The OBC and PBC results for r\ x coincide on the scale of the 
figure. Inset: A z = (SqSI /2 ) as a function of A > 1 (PBC) 
in comparison with analytic results. For A < 1.2 the length 
scale for saturation is larger than the system. 



fermions. For critical systems (see Fig. [|) we inves- 
tigated finite size effects on the estimates of the decay 
exponents rj x and r\ z (||) for < A < 1. For ordered 
chains with periodic boundary conditions DMRG results 
for (SqS^ 2 ) were in good agreement with the known 
long-range order parameter for A > 1 |24|] until the chain 
length became significantly shorter than the correlation 
(or saturation) length [^5| of the system (see inset of 
Fig.|l| ). These data suggest that chains of length L=80 
are sufficient to determine the long-range behavior of the 
correlation functions to better accuracy than could con- 
ceivably be obtained in the replica average. 

We then performed DMRG calculations for a large 
number of replicas each for various values of A and W. 
The number of replicas computed varied from 250 for 
small W, where replica expectation values fluctuate lit- 
tle, to more than 1000 for large W. We have not gath- 
ered replica averaged data for the XX-model by DMRG 
as the correlation functions for a small set of replicas 
showed perfect numerical agreement between DMRG and 
results from exact diagonalization for chain lengths of up 
to L = 160. 

Replica-averaged correlation functions in the critial 
and Ising regime are shown in Figure (^]) and (3) re- 
spectively. The data demonstrates qualitatively different 
behavior for small and large values of the disorder W. 
For A < 1 and W < 1 both x and z correlation func- 
tions decay algebraically with exponents that depend on 
A and W. Fitted to additional exponential components 
and finite offsets \ {SfSf) \ =A+\i-j\~ r i° exp(— -y|* — j j) 
the data show negligible inverse correlation lengths 7 and 
offsets A for small W and A. The decay of the x corre- 
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FIG. 2: 



Doubly logarithmic plots of replica averaged (a) 

{SgSr} and (b) (SqS*) for A = 0.5 and various disorder 
strengths W . z-correlation functions decay algebraically, the 
degree of correlation increases with disorder strength, while x- 
correlations decay algebraically for W < 1, but exponentially 
for W > 1. 



lation accelerates with growing disorder W whereas that 
of the z correlation decelerates. The values of the decay 
exponents r\ x and rj z as a function of W are shown in 
Figure (H)(a), indicating a continuous, but small change 
of both exponents from their values for ordered chains 
W = 0. 

For W > 1 the x-correlation functions decay expo- 
nentially in both the Ising and the critical regime. The 
inverse correlation length of the x correlation function 
is shown in Figure (|J)(b). The data is consistent with 
a crossover to exponential decay at W = 1 with signifi- 
cant finite size effects for W > 0.8, in particular in the 
Ising regime. In the Ising phase (A > 1) the z correla- 
tion functions continue to decay to a finite value for large 
separations (see Figure (||) (inset)). 

The nature of the transition at W — 1 is explained by 
a simple exactly solvable model, defined by the bimodal 
distribution 



p(X)=pS(X + l) + (l-p)S(X-l) 



(3) 



in the general Hamiltonian ([!]) with A < 1. The ground 
state spin correlations of this model are related to the 



known correlations |23 of the homogeneous chain (p = 0) 
by a simple gauge transformation p6| , |2?| . Consider a sin- 
gle replica, i.e. one configuration of A^s drawn from the 



3 




FIG. 3: Logarithmic plots of replica averaged correlation 
functions in the Ising regime (A = 1.5) for various disorder 
strengths W . x-correlations decay algebraically for W < 1 
and exponentially for W > 1, z-correlation functions satu- 
rate at a finite value that increases with the disorder strength 
(inset) 



distribution for an open chain. By a suitable prod- 
uct U of 7r rotations exp(2inSf ) about the local S z spin 
axes the disorder may be gauged away, i.e. H = UHW 
describes a configuration with A-j = 1. As the z spin 
components are invariant under U the z correlations of 
the disordered system are identical to those of the ho- 
mogeneous system. In contrast, a product of two x spin 
components aquires a string of random signs: 



i+r— 1 



(4) 



l=i 



The disorder average of simply yields the x correlation 
of the homogeneous system, multiplied by (1 — 2p) r . The 
decay of the z correlation thus remains algebraic, whereas 
the x correlation function is modified by an exponentially 
decaying factor, with a correlation length diverging with 
critical exponent equal to unity at the two quantum crit- 
ical points p c = 0, 1: £ x ~ 2|p-p \ ■ Applied to the prob- 
ability density used in the DMRG calculations, this ar- 
gument yields p = (W — l)/(2W) and = \nW for 
W > 1 (heavy solid line in Fig. 0(b)). For A = data 
for larger systems are available p6[ and the crossover 
from = to = lnT^ is visible more clearly. 

In the limiting cases that are accessible by alternate 
techniques, our results are in good agreement with both 
exact data for W — (Fig. Q) and numerical diagonal- 
ization results |2(| for long (N < 256) XX chains (see 
also |ll|). The latter also display clear deviations from 
the IRFP behavior predicted by the SDRG. For the z 
correlation an r -2 decay with the corresponding finite- 
size scaling behavior p8j remains perfectly intact from 
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FIG. 4: (a) Exponents n for the algebraic decay of the x and z 
correlation functions in the critical regime A < 1. In contrast 
to the predictions of the RSRG the exponents change contin- 
uously and slowly from their (finite system) values with no 
disorder. Circles, diamonds and squares denote A = 0.2, 0.5, 
and 0.8, respectively. Open symbols correspond to n z and 
filled symbols to n x ■ (b) Inverse correlation length for the de- 
cay of the x correlation functions for various A as a function 
of W. For W > 1 all data are consistent with the results for 
the exactly solvable model (heavy line). 



W = up to W = 2. The static z structure factor is 
linear in the wave vector q and independent of W. In 
contrast, the x correlation does not show finite-size scal- 
ing and the static x structure factor is neither linear nor 
VF-independent. The x correlation decays progressively 
faster as W grows. The data may be fitted to a power 
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law as long as W < 1 , but with an exponent significantly 
smaller than the value of 2 predicted for the RS phase. 
For W > 1 the decay is exponential, with an inverse cor- 
relation length proportional to the fraction of negative 
As (as in the exactly solvable model above). 

Combined, these results suggest a qualitative revision 
of the influence disorder is thought to have on quantum 
spin chains. No signs of attraction to the IRFP (with 
universal and isotropic algebraic decay of the spin corre- 
lation functions) predicted by the SDRG were observed 
in our study. Instead we observe a disorder-driven quan- 
tum phase transition at W = 1 for < A < 1, from 
a spin liquid phase with algebraically decaying correla- 
tions (with non-universal exponents) at W < 1 to a dif- 
ferent phase with exponential decay of the x correlations. 
These observations suggest to critically reexamine the ap- 
plicability of the SDRG to this system and to investigate 
the possible existence of a crossover length-scale beyond 
which the IRFP emerges as relevant ^9|. We hope that 
research on thermodynamic properties and spin corre- 
lations, in particular of the numerically much more ac- 
cessible XX chains, combined with the investigation of 
pre-asymptotic behavior in the the SDRG on the other 
hand will help resolve this issue. 

The mechanism for the destruction of the critical phase 
by disorder elucidated in this study may have interesting 
implications for experiments. For nonbounded disorder, 
where fluctuating signs of the couplings are present with 
varying probability for any W, one may anticipate the 
existence of a crossover length scale where algebraic de- 
cay crosses into exponential decay as a function of sys- 
tem size that may be observable by studying long, but 
finite chains. DMRG studies for Gaussian disorder, for 
which replica averages are much more difficult to con- 
verge, are presently under way to explore this scenario. 
The results of the present study are directly applicable 
to systems where a crystal field or easy plane generates a 
natural anisotropy. Extrapolating from both limits, they 
suggest the continued existence of critical behaviour for 
weak isotropic disorder in the isotropic model, a scenario 
that we will investigate more thoroughly in the future. 
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